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PREFACE 


The  Climatological  Probability  of  Cloud-Free  Line-of-Sight  (CPCFLOS)  program  was  developed  by  the 
United  States  Aii  Force  Environmental  Technical  Applications  Center's  Environmental  Simulation  and 
Techniques  Branch  (USAFETAC/SYT)  to  meet  the  many  requests  for  cloud-free  line-of-sight  (CFLOS) 
statistics  needed  in  the  testing  and  deployment  of  various  Department  of  Defense  optical  sensing  systems. 
The  CPCFLOS  program  merged  existing  capabilities  of  several  programs  (ECS3DCFL  and  DNYCPCF) 
into  one  and  corrected  some  inconsistenceies  found  in  these  individual  programs  under  project  number 
7081503. 

The  purpose  of  this  technical  note  is  to  describe  the  methodologies  and  algorithms  used  to  generate 
climatological  probabilities  of  CFLOS.  The  CPCFLOS  software  employs  Jiree  main  algorithms:  the 
Burger  Aerial  Algorithm  to  enhance  the  fidelity  of  observed  cloud  cover,  the  SRI  model  to  relate  fractional 
sky  cover  to  CPCFLOS,  and  the  Keilson-Ross  Algorithm  to  estimate  the  duration  probabilities  of  CFLOS. 
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1.  INTRODUCTION 


The  United  States  Air  Force  Environmental 
Technical  Applications  Center  (USAFETAC)  is 
frequently  tasked  to  provide  statistics  on  the 
probability  of  cloud-free  line-of-sight  (CFLOS)  for 
the  design  and  testing  of  the  optical  components  of 
various  Department  of  Defense  systems.  The 
climatological  probability  of  CFLOS  is  the 
percentage  of  time  a  cloud-free  line-of-sight  occurs 
at  a  particular  location  as  a  function  of  month, 
time-of-day,  and  viewing  angle.  However, 
conventional  weather  observations  only  contain 
data  on  cloud  cover,  not  on  the  distribution  of 
clouds  in  the  sky  dome.  As  a  result,  a  simulation 
model  is  required  to  relate  cloud  cover  climatology 
to  statistics  on  the  climatological  probability  of 
CFLOS.  CFLOS  calculations  can  be  required  for 
a  wide  variety  of  viewing  configurations.  This 
report  is  limited  to  one  class  of  viewing;  a  I  ne-of- 
sight  that  runs  from  a  point  on  the  earth's  -..urface 
to  a  point  in  space.  (For  discussions  of  other  types 
of  viewing  see  Harms,  1986;  V/arren,  1991a;  Rupp 
and  Warren,  1992). 


At  USAFETAC,  cloud  climatologies  are  derived 
from  two  main  sources.  The  first  is  the 
USAFETAC  surface  database,  which  is  a 
worldwide  collection  of  surface  weather 
observations.  The  second  data  source  is  the  Air 
Weather  Service  Real  Time  Nephanalysis 
(RTNEPH).  This  is  a  worldwide  cloud  analysis 
derived  from  surface  and  satellite  cloud  cover 
analyses  (Kiess  and  Cox,  1988).  From  either  of 
these  sources,  the  frequency  distributic.i  of  total 
cloud-cover  can  be  computed.  This  report 
describes  how  the  climatological  frequency 
distribution  of  cloud-cover  is  obtained  and  used  to 
compute  the  climatological  probability  of  CFLOS. 
The  probabilities  can  be  estimated  for  an 
instantaneous  point  in  time  or  for  a  specified  time 
window  (i.e.,  the  probability  of  a  continuous 
CFLOS  for  a  time  period  of  t  minutes).  The 
procedures  outlined  in  the  report  are  used  by  the 
USAFETAC  computer  program  CPCFLOS. 
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2.  CLOUD-COVER  FREQUENCY  DISTRIBUTION 


2.1  Introduction.  The  climatological  probability 
of  cloud-free  line-of-sight  (CPCFLOS)  is  derived 
from  the  climatological  cloud-cover  frequency 
distribution.  These  distributions  can  be  obtained 
from  surface  observations,  satellite  analyses,  or  a 
combination  of  these  two  sources.  None  of  these 
is  an  ideal  data  source.  S’  .race  observations  suffer 
several  problems;  e.g.,  observer  bias,  poor 
resolution  of  reporting  codes,  and  differences  in 
local  reporting  practices.  Satellite  analyses  have 


documented  deficiencies  as  well.  (For  a  review  of 
some  of  the  biases  in  the  RTNEPH  total  cloud- 
cover  analysis,  see  Lowther  et  ai,  1991).  More 
fundamentally,  surface  and  satellite  observations  of 
total  cloud-cover  are  actually  two  different 
physical  phenomena  (Snow,  1990).  This  point  is 
illustrated  in  Figure  1.  Although  no  conection  is 
made  for  this  fact  with  the  RTNEPH,  substantial 
errors  can  result  in  estimating  CPCFLOS  if  it  is 
not  taken  into  account. 


Figure  1.  Viewing  perspectives  of  total  cloud-cover.  S-  apparent  cloud  cover,  C-  cloud  cover, 
G  -  ground  cover,  N  -  sky  cover.  C  is  the  view  angle  (from  Snow,  1990). 


These  problems  must  be  all  taken  into  account  to 
properly  compute  CPCFLOS  probabilities.  As  a 
result,  the  computation  of  the  cloud-cover 
frequency  distribution  is  not  as  straightforward  as 
merely  binning  cloud-cover  reports  over  a  given 
period  of  record.  A  statistical  technique  must  be 
used  to  convert  the  low  fidelity  cloud-cover  reports 
to  a  higher  fidelity  distribution.  Separate  binning 


must  be  made  of  surface-based  and  satellite-based 
cloud-cover  reports.  This  is  very  important  when 
processing  RTNEPH  cloud-cover  reports  since  the 
RTNEPH  is  a  blend  of  these  two  sources  of  data. 
At  some  locations,  the  climatology  will  be  nearly 
entirely  surface  or  satellite-based;  at  other 
locations  it  will  be  a  mixture  of  these  two  sources. 
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22  Surface  Reports  of  Total  Cloud  Cover.  obstruction  to  the  cloud  layer  (fog,  haze,  and 

Surface  observations  are  reported  with  either  the  snow  are  probably  the  most  frequently  occurring 

Airways,  METAR,  or  synoptic  code.  Each  of  obscurants).  A  total  obscuration  means  that  the 

these  codes  uses  a  different  encoding  procedure  to  entire  sky  dome  is  obstructed.  If  any  portion  of  a 

report  the  observations.  cloud-layer  o-  the  clear  sky  is  visible,  then  the 

obstruction  is  partial.  The  observer  then  carries  a 
2.2.1  Airways  Code.  With  tne  Airways  code,  the  report  in  the  remarks  section  of  the  observation 

observer  reports  cloud  cover  as  one  of  four  stating  the  fraction  (in  tenths)  of  the  sky  dome  that 

categories:  clear  (CLR),  scattered  (SCT),  broken  is  obscured  and  the  nature  of  the  obscurant  (e.g., 

(BKN),  or  overcast  (OVC).  The  observer  reports  the  remark  'T8"  means  that  fog  obscures  8/10  of 

the  total  cloud  cover  and  base  of  each  layer  of  the  sky  dome).  Cloud  distributions  for  airways 

cloud,  beginning  with  the  lowest  layer.  The  reporting  stations  are  initially  obtained  by 

reported  total  cloud  cover  reported  is  cumulative.  computing  the  relative  frequency  of  each  of  the 

Thus  the  category  reported  for  the  last  layer  (the  four  cloud  classes.  Generally,  only  a  small 

layer  with  the  highest  observed  base),  is  also  the  percentage  of  observations  are  reported  as  either 

total  cloud-cover.  Unfortunately,  there  is  not  a  totally  obscured  or  partially  obscured.  Totally 

common  standard  defuiition  of  these  four  obscured  conditions  can  be  binned  with  the  OVC 

categories.  Table  1  compares  the  Air  Force/Navy  category.  Partial  obscuration  of  the  sky  dome  can 

and  National  Weather  Service/Federal  Aviation  be  as  small  as  one-tenth  or  as  large  as  nine-tenths. 

Administration  definitions  of  these  categories  We  have  arbitrarily  selected  to  bin  these  partial 

(FMH-1,  1988).  In  addition  to  these  four  cloud-  obscurations  as  BKN.  However,  the  fraction  of 

cover  categories,  two  other  classes  may  be  observations  reporting  a  partial  obscuration  is 

reported:  partially  obscured  and  totally  obscured.  small,  and  this  assumption  will  have  little  (if  any) 

An  obscuration  is  a  surface-based  visual  effect  on  the  subsequent  frequency  distribution. 


TABLE  1.  Comparison  of  the  definitions  of  tho  Ainwa/s  cljud-cover  categories  as  used  in  the 
United  States  (FMH-I,  1988). 


Category 

National  Weather  Service/ 

Federal  Aviation  Administration 

Air  Force/Navy 

CLR 

Cloudless  or  sky-cover  less  than  1/10. 
(round  to  the  nearest  tenth) 

Cloudless. 

SCT 

Sky-cover  between  1/10  and  5/10 
inclusive. 

Sky-cover  between  a  trace 
and  5/10,  inclusive. 

BKN 

6/10  through  9/10,  inc'usive. 

Greater  than  or  equal  to 

6/10,  but  not  completely 
overcast. 

OVC 

Greater  than  9/10,  or  completely 
overcast. 

Completely  overcast. 

2.2.2  Synoptic  Code.  The  synoptic  code  has  a 
field  reserved  for  reporting  the  total  cloud-cover  in 
eighths  (a  single  number  between  0  and  8  is  used). 
As  with  the  Airways  code,  two  additional  values 
are  used  for  obscured  conditions:  9  for  partially 
obscured  and  10  for  totally  obscured.  We 
converted  the  totally  obscured  reports  to  8 
(overcast)  and  partially  obscured  to  4  (four- 
eighths). 


cover  reports.  Separate  frequency  distributions 
should  then  be  computed  for  surface-  and  satellite- 
based  reports.  While  the  type  of  surface  report  is 
not  included,  it  is  fairly  easy  to  tell  from 
examination  of  the  frequency  distribution  if  the 
surface  observations  are  primarily  Airways  (large 
number  of  counts  at  C5fc,  25%,  85%,  and  100% 
total  cloud-rover),  or  METAR/synoptic  (counts  at 
0.  12%,  25%,  38%,  etc.) 


2.Z3  METAR  Code.  METAR  code  is  the  most 
difficult  to  work  with  in  computing  the  total  cloud- 
cover  distribution.  This  is  due  to  the  fact  there  is 
no  provision  in  the  code  for  reponing  the  total 
cloud-cover.  Cloud  cover  is  reported  only  in 
layers,  but  these  layered  amounts  are  not 
cumulative  as  in  the  Airways  code.  The  observer 
reports  the  fraction  of  cloud  (in  eighths)  that  is 
observed  in  each  layer.  Merely  adding  the  layer 
coverage  does  not  provide  the  total  cloud  cover, 
since  the  observer  may  report  overlapping  layers. 
An  algorithm  has  been  developed  by  USAFETAC 
to  process  these  layered  amounts  to  obtain  a  total 
cloud-cover;  it  is  described  in  the  Appendix. 
Values  of  partial  and  total  obscuration  are 
processed  the  same  as  with  the  synoptic  code. 

2-3  RTNEPH  Total  Cloud  Cover.  The 
RTNEPH  total  cloud  cover  is  recorded  to  the 
nearest  whole  percent.  There  are  two  primary 
sources  of  these  reports:  surface  observations  and 
satellite  analyses.  If  a  nearby  surface  observation 
is  available  for  the  RTNEPH  point  of  interest,  this 
observation  is  used  with  no  additional  processing. 
The  corresponding  percentages  for  the  Airways 
code  categories  of  clear,  scattered,  broken,  and 
overcast  are  0%,  25%,  85%,  and  100%  Reports 
for  the  synoptic  code  are  converted  to  the 
corresponding  percentage.  RlTfEPH  uses  an 
algorithm  to  determine  the  total  cloud  cover  from 
METAR  reports.  If  no  surface  observations  are 
available,  then  a  satellite  analysis  is  used  and 
reported  in  whole  percent  (for  details,  see  Kiess 
and  Cox,  1988).  The  RTNEPH-based  total  cloud 
cover  distribution  will  then  be  a  mix  of  these 
procedures.  There  is  a  flag  in  the  RTNEPH 
database  indicating  the  source  of  the  total  cloud 


2.4  Burger  Aerial  Algorithm.  Because  of  the 
discrete  nature  of  surface  cloud-cover  reports,  we 
make  use  of  the  Burger  Aerial  Algorithm  (BAA) 
to  convert  the  observed  distribution  to  one  with 
higher  fidelity.  The  BAA  (Burger,  1985)  is 
actually  a  continuous  distribution  of  cloud-cover, 
but  for  practical  reasons,  we  use  it  to  determine 
the  probability  of  each  of  the  cloud-cover 
frequency  intervals  listed  in  Table  2.  (For  a 
comprehensive  review  of  this  algorithm  and  a 
discussion  of  its  accuracy,  see  Henderson-Sellers 
and  McGuffie,  1991).  The  computational  details 
of  converting  the  Airways,  METAR,  or  synoptic 
code  to  these  intervals  using  the  BAA  are 
described  in  the  Appendix.  The  categories  in 
Table  2  are  used  to  compute  the  climatological 
probability  of  CFLOS,  as  described  in  the  next 
chapter. 

TABLE  2.  Cloud-cover  intervals  used  for 
CPCFLOS  calculations. 


Category 

Range 

Number 

{%) 

2  6'' 

15" 

•J  ! 

-1  26-' 

30': 

5  36-~ 

-  45-7 

5  46' 

7  56^ 

~  -JJ  z 

r'  r.  ^ 

9  76" 

*■  J  56' 

^ 

-  55^ 

_ 

:  1  ^6^ 

-  !0C^ 
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3.  CALCULATION  OF  THE  CLOUD-FREE  LINE-OF-SIGHT  PROBABILITIES 


3.1  SRI  Model.  Lund  and  Shanklin  (1973) 
developed  a  technique  for  relating  fractional  sky- 
cover  to  the  climatological  probability  of  cloud- 
free  line-of-sight.  This  technique  was  derived  from 
a  collection  of  whole  sky  photos  (Lund  and 
Grantham,  1980)  taken  at  Columbia  MO.  Marick 
et  al.  (1979)  identified  shortcomings  of  the  Lund 
and  Shanklin  model  and  proposed  a  revised  model 
that  is  partly  derived  from  relationships  in  the 
cloud  viewing  geometry.  This  revised  model  is 
referred  to  as  the  Stanford  Research  Institute  (SRI) 
model.  Malick  et  al.  (1979)  describes  this  model 
in  detail.  In  the  SRI  model,  the  probability  of  a 
cloud-free  line-of-sight  when  looking  straight  up, 
is  given  by: 

=  1  -  (1) 

4 

where  S  is  the  fractional  sky-cover  as  viewed  from 
the  ground.  For  other  viewing  angles,  the 
following  approach  is  used.  First  we  define  the 
zenith  angle,  6,  as  the  angle  the  desired  line-of- 
sight  makes  with  the  vertical.  The  climatological 
probability  of  CFLOS,  P(Q^,  given  a  zenith  angle 
of  0  and  a  fractional  sky-cover  5,  is  obtained 
using: 

F\Q,S)  =  P^5)<’  *  (2) 

where  b  is  the  average  height-to-width  ratio  of 
clouds  and  is  assumed  to  follow  the  relation: 

b  =  0.55  -  I  (3) 

To  obtain  the  climatological  probability  of  CFLOS 
for  viewing  angle  0  [FjC®)]  the  CFLOS  probability 
for  a  given  sky -cover  value  S  [F(0,S)]  is  weighted 
by  the  relative  frequency,  f(S),  of  the  sky-cover 
value: 


pjie)  =  fr^s)p(d.S)dS  W 

0 

In  practice,  a  finite  approximation  to  Equation  4  is 
used: 

M 

where  /  is  the  relative  frequency  of  cloud-cover 
interval  i  and  P,(6)  is  the  corresponding  CFLOS 
probability  for  this  interval.  The  intervals 
typically  used  are  shown  in  Table  2. 

A  comparison  between  observed  cloud-free  line- 
of-sight  from  GOES  imagery  with  computed 
values  from  the  SRI  model  are  shown  in  Figure  1. 
The  results  show  very  good  agreement.  Warren 
(1991b)  also  presents  a  favorable  comparison 
between  predicted  (using  the  SRI  model)  and 
observed  climatological  probabilities  of  sunny  line- 
of-sight  in  Hawaii. 


Figure  2.  Comparison  between  predicted 
CPCFLOS  values  (from  the  SRI  model)  and 
values  derived  from  GOES  imagery  at  three 
locations  (Results  courtesy  of  Al  Boehm, 
Hughes  STX  Corporation). 
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3^  Space-Based  Observations  of  Clouds.  The 
SRI  model  assumes  the  fractional  cloud-cover  is 
deterrr.uied  by  a  surface  observer.  For  locations 
with  no  surface  observations  available,  cloud 
climatologies  must  be  derived  from  satellite-based 
instruments.  Such  an  approach  is  used  by  the 
USAF  Real  Time  Nephanalysis  (RTNEPH). 
However,  surface  and  satellite  observations  of 
cloud-cover  are  distinct  phenomena  (Snow,  1990). 
Surface-based  observations  of  fractional  cloud- 
cover  are  often  referred  to  as  "sky-cover,*  while 
space-based  observations  are  referred  to  as  "earth- 
cover."  The  use  of  the  SRI  model  with  a  cloud 
climatology  derived  from  satellite  data  may 
produce  significant  errors.  (A  recent  unpublished 
USAFETAC/SYT  study  found  errors  as  high  as 
26%,  with  typical  errors  of  10-15%). 

If  we  assume  that  the  space-based  observation  is 
obtained  by  looking  straight  down  at  the  point,  the 
corresponding  earth-cover  for  a  given  value  of  sky- 
cover  can  be  obtained  from: 

O  =  1  -  (6) 


where  Q  is  the  fractional  earth-cover,  and  P„  is  the 
probability  of  a  surface-based  cloud-free  line-of- 
sight  for  vertical  viewing.  Combining  Equation  4 
with  Equation  2,  the  corresponding  sky-cover,  (S), 
for  a  given  value  of  earth-cover,  (^,)  is: 


S  = 


2 

6 


(7) 


(NOTE:  This  equation  is  valid  for  climatological 
use  only;  in  any  given  cloud  scene,  the 
relationship  between  S  and  Q  must  be  defined  in 
terms  of  a  statistical  distribution).  For  a  given 
earth-cover,  Q,  using  Equation  7  we  can  obtain  the 
corresponding  sky-cover,  5,  and  then  compute 
climatological  CFLOS  probabilities  with  Equations 
1  through  3.  Table  3  lists  CFLOS  probabilities  for 
surface-based  observations.  Table  4  is  the 
corresponding  table  for  satellite-based  cloud-cover. 

33  Combining  Sources.  For  RTNEPH-based 
calculations  of  CPCFLOS,  two  frequency 
distributions  may  be  available  for  a  point  in 
question:  one  based  on  surface  reports  and  one 
based  on  analyses  of  satellite  observations. 
Separate  CFLOS  probabilities  are  obtained  from 
each  source.  The  final  computed  value  is  just  a 
weighted  average  of  the  two  sources: 

CPCFLOS  =  (8) 

P,  ♦  Pj 

where  F,  and  F^  are  the  CPCFLOS  probabilities 
determined  front  surface  data  and  satellite  dau 
respectively,  and  Nj  and  Nj  are  the  corresponding 
number  of  observations  used  for  each  compulation. 
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TABLE  3.  CFLOS  probabilities  as  a  function  of  viewing  angle  for  each  of  the  eleven  cloud* 
coverintervals  defined  in  Table  1 .  {Surface-based  observations) 


Deg  Cloud-Cover  Category 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 

0.99 

0.97 

0.92 

0.86 

0.78 

0.69 

0.58 

0.46 

0.32 

0.17 

0.04 

10 

0.99 

0.96 

0.91 

0.85 

0.77 

0.67 

0.57 

0.45 

0.31 

0.16 

0.04 

20 

0.99 

0.96 

0.91 

0.84 

0.76 

0.66 

0.55 

0.43 

0.30 

0.16 

0.04 

30 

0.99 

0.96 

0.90 

0.83 

0.74 

0.64 

0.54 

0.42 

0.29 

0.15 

0.04 

40 

0.99 

0.95 

0.89 

0.81 

0.73 

0.63 

0.52 

0.40 

0.28 

0.14 

0.04 

50 

0.99 

0.95 

0.88 

0.80 

0.70 

0.60 

0.49 

0.38 

0.26 

0.14 

0.03 

60 

0.99 

0.94 

0.86 

0.77 

0.67 

0.57 

0.46 

0.35 

0.24 

0.12 

0.03 

70 

0.98 

0.92 

0.83 

0.72 

0.61 

0.50 

0.40 

0.30 

0.20 

0.10 

0.03 

80 

0.97 

0.88 

0.74 

0.61 

0.48 

0.36 

0.27 

0.19 

0.12 

0.06 

0.01 

TABLE  4.  CFLOS  probabilities  as  a  function  of  viewing  angle  for  each  of  the  eleven  cloud-cover 
inten/als  defined  in  Table  1.  (Satellite-based  observations) 


Deg  Cloud-Cover  Category 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 

0.98 

0.90 

0.80 

0.70 

0.60 

0.50 

C40 

0.30 

0.20 

0.10 

0.02 

10 

0.97 

0.89 

0.79 

0.69 

0.59 

0.49 

0.39 

0.29 

0.19 

0.10 

0.02 

20 

0.97 

0.89 

0.78 

0.67 

0.57 

0.47 

0.38 

0.28 

0.19 

0.10 

0.02 

30 

0.97 

0.88 

0.76 

0.66 

0.56 

0.46 

0.36 

0.27 

0.18 

0.09 

0.02 

40 

0.96 

0.87 

0.75 

0.64 

0.54 

0.44 

0.35 

0.26 

0.17 

0.09 

0.02 

50 

0.96 

0.85 

0.73 

0.61 

0.51 

0.42 

0.33 

0.24 

0.16 

0.08 

0.02 

60 

0.95 

0.83 

0.70 

0.58 

0.48 

0.39 

0.30 

0.22 

0.15 

0.07 

0.02 

70 

0.94 

0.79 

0.64 

0.52 

0.42 

0.33 

0.26 

0.19 

0.12 

0.06 

0.01 

80 

0.90 

0.69 

0.51 

0.38 

0.28 

0.21 

0.16 

0.14 

0.07 

0.04 

0,01 

4.  PROBABILITY  OF  CFLOS  DURATION 


4.1  Keilson-Ross  Procedure.  The  technique 
described  in  Chapter  3  provides  the  probability  of 
CFLOS  for  an  instant  in  time.  For  some  systems, 
a  CFLOS  window  for  successful  viewing  may  be 
required.  A  technique  for  estimating  duration 
probabilities  of  meteorological  conditions  is 
described  by  Gringorten  (1982)  and  summarized 
here.  This  technique  is  known  as  the  "Keilson- 
Ross  algorithm"  (Ross,  1980).  Hering  (1989)  used 
this  approach  for  computing  the  climatological 
probability  of  CFLOS  duration. 

If  the  probability  of  an  event  depends  solely  upon 
the  outcome  of  the  preceding  trial,  the  sequence  of 
events  is  called  a  Markov  process,  after  the 
Russian  mathematician  A. A.  Markov  (Frieden, 
1983).  Weather  elements  are  often  modeled  as 
Markov  processes.  The  Omstein-Uhlenbeck  (O-U) 
process  is  a  class  of  Markov  processes  in  which 
the  correlation  coefficient  (p)  over  the  time 
interval  (At)  between  the  equivalent  normal 
deviates  (ENDs)  of  the  two  endpoints  (y^  yt.4,); 


where  t  is  a  parameter  known  as  the  "relaxation 
time."  Hering  (1989;  1990)  recommends  a  x- 
value  of  30  minutes  to  calculate  the  correlation  of 
CLOUD/NO  CLOUD  as  a  function  of  time  for  a 
fixed  point  in  the  sky  (p).  (For  a  discussion  of 
equivalent  normal  deviates,  see  section  A.2  in  the 
Appendix.) 

Given  the  value  of  a  weather  element  at  time  t 
with  a  conesponding  END  of  an  END  at  time 
t  -r  At  can  be  obtained  stochastically  using 
(Whiton  and  Berecek,  1982): 

y  -  py,  *  T1V1  -  P* 
where  rj  is  a  random  normal  number. 


The  cumulative  probability  that  a  weather  element 
jc  is  less  than  or  equal  to  a  threshold  value  (x^)  is 
the  same  probability  that  its  ENT)  (y)  is  less  than 
or  equal  to  the  END  of  the  threshold  value  (y^). 
This  probability  can  thus  be  obtained  from 
integrating  the  normal  distribution  curve: 

P{x  i  =  P[y  i.  yi 


We  can  use  the  notation  F^/x^)  to  represent  the 
cumulative  probability  of  a  weather  element 
having  a  value  less  than  or  equal  to  a  threshold 
(Xj.)  throughout  a  time  duration  Af.  This 

probability  cctresponds  to  the  probability  that  the 
END  of  the  v.  ither  element  will  have  a  value  less 
than  or  equal  to  yy  END  of  for  a  relative 
time  period  0  (this  is  a  dimensionless  value--for 
the  O-U  process  this  corresponds  to  Ar/r): 

FM  -  ^a(F7)  (^2) 

Equations  10  and  1 1  can  be  combined  to  obtain  a 
functional  relationship  for  F^Cy^.).  For  the  special 
case  when  yy  =  0  (i.e.,  the  case  when  Xj  is  the 
median  value): 

Fe(0)  =  -sin-U©-^  03) 

Unfortunately  this  simple  result  is  limited  only  to 
the  specific  case  when  the  threshold  is  the  median 
value.  An  analytical  expression  for  the  general 
solution  ha<;  not  been  found.  A  complete 
description  of  the  mathematics  involved  in 
obtaining  an  iterative  solution  is  described  in 
Gringorten  (1982).  Ross  (1980)  fitted  this  solution 
to  a  series  of  polynomials  obtained  from  the  use  of 
a  cubic  spline  algorithm.  Gringorten  (1982) 
provides  the  FORTRAN  code  to  implement  this 
algorithm. 
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A2  Example.  Suppose  at  a  given  location  and 
time,  the  climatological  probability  of  CFLOS  (at 
an  instant  in  time)  tor  a  zenith  angle  of  40®  is 
49%.  We  can  use  the  Keilson-Ross  procedure  to 
compute  the  probability  of  CFLOS  for  time  periods 
of  t  minutes.  To  plot  the  CPCFLOS  as  a  function 
of  the  time  window,  we  can  first  compute  the 
values  at  selected  points,  such  as  1,  2,  5,  10,  20, 
30, 45,  and  60  minutes.  The  correlations  obtained 
from  the  O-U  assumption  are  listed  in 


TABLE  5.  Correlation  between  ENDs  of 
cloLid/no  cloud  values  (p)  and  corresponding 
CPCFLOS  as  function  of  time  window. 


Time 

9 

Probability 

0 

1.00 

0.49 

1 

0.97 

0.46 

2 

0.94 

0.45 

5 

0.85 

0.42 

10 

0.72 

0.39 

20 

0.51 

0.35 

30 

0.37 

0.33 

45 

0.22 

0.29 

60 

0.14 

0.27 

Table  5.  Also  in  Table  5  are  the  probabilities  of 
CFLOS  for  the  specified  duration  obtained  from 
the  Keilson-Ross  algorithm.  A  plot  of  the  results 
is  shown  in  Figure  3.  Results  of  a  validation 
study  for  this  technique  in  CFLOS  calculations  are 
presented  in  Hering  (1989;  1990).  This 

assumption  that  cloud-free  line-of-sight  is  assumed 
to  be  a  Markov  process  is  clearly  the  limitation  of 
this  technique.  Other  models  have  been  suggested 
(McCabe,  1968;  Boehm,  19926). 


sight  as  a  function  of  duration  window. 
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5.  SAMPLE  CALCULATION 


To  illustrate  how  the  climatological  probability  of 
CFLOS  is  obtained,  consider  the  following 
example.  Table  6  gives  an  example  with  two 
sample  cloud-cover  distributions.  From  this  data, 
the  parameters  of  the  Burger  Aerial  Algorithm  are 
computed  using  the  method  described  in  Appendix 
A.  For  this  example,  the  mean  sky  cover  values 
are  0.733  and  0.490,  with  sky  dome  scale  distances 
of  6.551  km  and  1.604  km  for  Points  A  and  B 
respectively.  The  relative  frequency  distribution  of 


TABLE  6.  Two  sample  cloud-cover 
frequency  distributions. 


Number  of 

Relative 

Observations 

Freauencv 

Point  A 

CLR  346 

0.208 

SCT  101 

0.061 

BKN  98 

0.058 

OVC  1120 

0.673 

Point  B 

CLR  414 

0.263 

SCT  376 

0.248 

BKN  412 

0.258 

OVC  364 

0.231 

cloud-cover,  using  the  Burger  intervals,  are 
presented  in  Table  7.  From  this  distribution,  the 
climatological  probability  of  CFLOS  can  be 
obtained  using  the  data  in  Table  3.  These  results 
are  presented  in  Table  8.  The  first  column  of 
results  is  for  a  CFLOS  window  of  0  (instantaneous 
probability),  the  second  column  is  for  a  window  of 
10  min.  These  latter  numbers  were  obtained  by 
first  using  the  Keilson-Ross  algorithm  to  adjust  the 
probabilities  of  CFLOS. 


TABLE  7.  Relative  cloud-cover  frequency 
distribution  for  1 1  cloud-cover  categories. 


Cloud-Cover 

Relative  Freauencv 

A 

B 

<  0.05 

0.21 

0.27 

0.06  -  0.15 

0.02 

0.08 

0.16  -  0.25 

0.01 

0.05 

0.26  -  0.35 

0.01 

0.05 

0.36  -  0.45 

0.01 

0.04 

0.46  -  0.55 

0.01 

0.04 

0.56  -  0.65 

0.01 

0.04 

0.66  -  0.75 

0.01 

0.05 

0.76  -  0.85 

0.01 

0.05 

0.86  -  0.95 

0.02 

0.08 

0.96  -  1.00 

0.68 

0.25 
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TABLE  8.  Climatological  Probability  of  CFLOS  for  each 
of  the  two  sample  datasets. 


Zenith 

10-minute 

Anale 

Instantaneous 

Window 

Point  A 

0 

31 

26 

10 

31 

26 

20 

30 

26 

30 

30 

25 

40 

29 

25 

50 

29 

25 

60 

29 

24 

70 

28 

24 

80 

26 

22 

Point  B 

0 

58 

48 

10 

58 

48 

20 

57 

47 

30 

57 

47 

40 

56 

46 

50 

55 

45 

60 

54 

44 

70 

52 

42 

80 

47 

37 
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6.  PROGRAM  STRUCTURE 


The  techniques  described  in  Chapters  3  and  4  are 
used  by  the  USA.reTAC  CPCFLOS  program.  The 
program  ingests  cloud  data  from  several  possible 
sources  (surface  observations;  the  old  3-D 
nephanalysis,  or  SDNEPH;  and  the  real-tim' 
nephanalysis,  or  RTNEPH)  and  derives  the 


climatological  probability  of  cloud-free  line-of- 
sight.  A  diagram  of  the  program  structure  is 
depicted  in  Figure  4.  The  first  row  of  subroutines 
snown  in  this  figure  lists  the  various  readers  used 
by  CPCFLOS,  depending  on  the  type  of  data 
selected. 


CPCFLOS  uses  one  of  four  readers  to  process  sky- 
cover  data.  Subroutine  DB2RETV  retrieves  data 
stored  in  DB2  tables  (DB2  is  the  relational 
database  management  system  on  the  USAFETAC 
mainframe.  SFCRETV  is  used  for  DATSAV2 
surface  data  (DATSAV2  is  a  USAFETAC- 
developed  code  format  for  archiving  surface 
observations;  for  details  see  USAFETAC,  1986), 
RTRETV  is  used  for  RTNEPH  data,  and 
D3RETV  for  3DNEPH,  an  earlier  version  of 
RTNEPH;  for  details  see  Cobum,  1970).  After 
each  observation  is  read  in,  the  data  is  parsed  to 
subroutine  BINCNT,  wv  ch  keep;  track  of  the 
number  of  occurrences  of  each  sk'  .o  .'er  category, 
by  month,  and  hour.  RTNEPH  data  is  also 
stratified  by  wiiether  the  cloud-cover  was  derived 
from  a  surface  report  or  a  satellite  analysis. 


The  RTNEPH  became  operational  on  1  June  1984. 
The  drawback  with  3DNEPH,  its  predecessor,  is 
that  no  data  source  flag  is  available;  users  cannot 
be  certain  if  an  observation  is  surface-  or  satellite- 
based.  CPCFLOS  includes  a  flag,  set  by  the  user, 
to  designate  if  there  is  a  surface  observing  station 
within  a  50-NM  radius  of  the  desired  point.  If 
such  a  station  exists,  CPCFLOS  assumes  that  the 
ZJj  lEPH  cloud-cover  reports  are  surface-based;  if 
not,  it  assumes  they  are  satellite-based.  Because 
of  this  serious  limitation,  the  use  of  3DNEPH  to 
compute  CFLOS  probabilities  should  be  avoided 
whenever  possible. 

Once  all  the  observations  have  been  read  in,  the 
subroutine  BINFRQ  is  called.  This  routine 
computes  the  relative  frequency  of  observed  cloud- 
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cover.  Categories  are  in  eighths  of  tota<  cloud 
cover.  (For  Airways  reports,  clear  skies  are 
assigned  a  value  of  0/8;  scattered,  2/8;  broken,  7/8; 
and  over-ast  8/8). 

Next,  control  is  passed  to  the  subroutine 
BAAPRM,  which  computes  the  parameters  of  the 
Burger  Aerial  Algorithm.  Two  techniques  are 
used:  one  for  Airways  data  and  one  for  data  in 
eighths.  An  algorithm  is  used  to  examine  the 
relative  frequency  distribution  passed  from 
BINFRQ.  If  the  total  number  of  observations  in 
2/8  exceeds  the  combined  total  in  1/8, 3/8,  and  4/8, 
and  if  the  total  number  in  the  category  7/8  exceeds 
the  combined  total  in  S/8  and  6/8,  the  distribution 
is  classified  as  Airways.  This  classification  may 
occur  with  both  surface-  and  NEPH-based  data. 

In  the  following  step,  the  relative  frequency 
distribution  of  cloud  cover  for  the  11  intervals 
presented  in  Table  2  is  computed  by  the  subroutine 
CLDFRQ,  which  calls  the  subroutine  "BURGER." 
When  passed  the  mean  sky  cover,  scale  distance, 
and  a  value  of  cloud-cover,  BURGER  returns  the 
probability  of  obtaining  a  cloud-cover  less  than  or 
equal  to  the  paissed  value. 


Next,  the  subroutine  PCFLOS  computes  the 
climatological  probability  of  CFLOS  given  the 
relative  frequency  distribution  computed  by 
CLDFRQ.  PCFLOS  begins  by  computing  the 
CFLOS  probabilities  for  each  cloud -cover  interval 
and  viewing  angle,  first  for  surface-based,  then  for 
satellite-based  data.  Next,  if  a  time  window  other 
than  zero  is  requested,  each  probability  is  passed 
to  the  function  DURATN,  which  adjnrts  the 
probability  for  the  time  window  using  the  Keilson- 
Ross  procedure.  This  procedure  is  implemented 
by  the  subroutine  RKEROS,  which  is  an  empirical 
fit  to  the  analytical  solution  using  cubic  splines. 
The  subroutine  BSPERVS  is  needed  to  solve  the 
cubic  spline  equations.  PCFLOS  returns  separate 
CFLOS  probability  estimates:  one  for  surface- 
based  data,  the  other  for  satellite-based  data. 

The  final  routine  to  be  called  is  MERGER,  which 
combines  the  two  CFLOS  probabilities  into  one, 
using  a  weighted  average.  The  individual 
probabilities  are  weighted  by  the  number  of 
observations  used  to  derive  it.  This  is  the  final 
value  of  CFLOS  probability,  which  then  appears 
on  the  output.  The  program  output  is  created  by 
subroutine  CFLOUT. 
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APPENDIX 


BURGER  AERIAL  ALGORITHM 


A.1  Introduction.  Surface  observations  of  total 
cloud-cover  are  reported  with  low  resolution.  In 
the  U.S.,  where  most  observations  are  recorded  in 
the  Airways  code,  cloud-cover  is  reported  as  either 
clear,  scattered,  broken,  or  overcast.  Selected 
locations  in  the  U.S.  (and  in  most  other  countries 
as  well)  record  an  observation  in  synoptic  code 
every  3  hours.  In  synoptic  units  are  in  eighths  of 
total  cloud-coverage.  At  most  other  locations 
throughout  the  world,  the  usual  code  for  weath  r 
observations  is  METAR.  Only  layer  amounts  are 
encoded  in  METAR  (in  eighths);  there  is  no  report 
for  total  cloud-cover. 

The  limitations  of  using  these  codes  to  infer  the 
distribution  of  cloud-cover  is  reviewed  by 
Grantham  and  Boehm  (1986).  To  overcome  the 
low  fidelity  of  the  observations,  Burger  (1985) 
devised  a  model  for  specifying  the  cloud-cover 
distribution.  We  refer  to  this  model  as  the  Burger 
Aerial  Algorithm  (BAA).  The  model  is  a 
numerical  fit  to  the  results  of  multiple  simulations, 
as  described  by  Burger  and  Gtingorten  (1983).  It 
converts  the  discrete  observations  present  in  all 
three  codes  to  a  continuous  distribution.  This 
algorithm  was  favorably  reviewed  by  Henderson- 
Sellers  and  McGuffie  (1991).  Other,  somewhat 
more  complex,  algorithms  have  been  proposed 
(Boehm,  1992a). 

The  BAA  has  two  fundamental  variables;  mean 
sky  cover  and  sky-dome  scale  distance.  Frequency 
distribution  for  various  values  of  these  variables 
are  shown  in  Figure  A-1,  opposite.  A  large  scale 
distance  indicates  cloud-cover  distributions  that  are 


either  predominantly  clear  or  cloudy.  In  these 
instances,  the  mean  is  the  least  likely  value  to 
occur.  A  small  scale  distance,  on  the  other  hand, 
signifies  a  distribution  centered  about  the  mean. 

A.2  Equivalent  Normal  Deviates.  To 
understand  the  fundamentals  of  the  Burger  Aerial 
Algorithm,  a  basic  understanding  of  transforming 
variables  to  a  Gaussian  domain  (a  process 
sometime  referred  to  as  transnormalization)  is 
required.  For  details,  see  Boehm  (1976)  or 
Whiton  and  Berecek  (1982).  A  normal  deviate  is 
a  variable  that  follows  a  normal  distribution  with 
a  mean  of  zero  and  a  standard  deviation  of  unity. 
A  value  of  any  other  variable  will  have  a  one-to- 
one  correspondence  with  this  normal  deviate 
through  the  cumulative  frequency  distribution. 
The  corresponding  value  of  the  normal  deviate 
(the  one  with  the  same  cumulative  probability  as 
the  value  of  the  meteorological  variable  of 
interest)  is  referred  to  as  its  equivalent  normal 
deviate  (END).  For  example,  assume  the  observed 
frequency  of  total  cloud-cover  of  three-eighths  or 
less  at  a  given  location  is  35%.  The 
corresponding  value  of  the  normal  deviate  with  a 
cumulative  probability  of  35%  is  -0.385.  In  this 
case,  the  END  of  three-eighths  cloud-cover  is 
-0.385.  This  is  an  example  of  transforming  to  the 
Gaussian  domain.  The  advantage  of  this 
transformation  is  that  it  permits  the  use  of 
statistical  operations,  which  require  an  assumption 
of  normality.  Most  meteorological  variables, 
especially  the  cloud-cover,  are  especially  non¬ 
normal. 
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Rgure  A-l.  Relative  frequency  distribution  of  doud-cover  for  various  values  of  Itie  rrean  doud-cover  and  sicy  dome 
scale  distance  (Burger,  1985). 


15 


A3  Model  Algorithm 

A.3.1  Airways  Code.  Cloud-cover  reports  in 
Airways  code  are  either  clear,  scattered,  broken,  or 
overcast.  If  we  define  the  relative  frequency  of 
each  category  as  CLR,  SCT,  BKN,  and  OVC,  the 
BAA  variables  are  computed  as  follows:  First,  the 
ENDs  of  CLR,  SCT,  and  BKN  are  computed 
(Boehm,  1976)  and  referred  to  as  y„,  and  y,^ 
Differences  between  these  values  are  then 
determined: 

(Ay)i  -  y»r-  y* 

=  /*»  -  /dr 

(Alls  =  /«„  -  /^ 

Next  we  consider  three  cases.  Case  1  will  be  valid 
for  most  locations.  Case  2,  by  exception,  is  when 
OVC  is  near  zero.  Case  3,  also  by  exception,  is 
when  CLR  is  near  zero.  The  scale  distance 
parameter,  (z),  is  computed  using: 

z  *  a  *  b  ln(A/  ♦  d)  05) 


TABLE  A-1.  Coefficients  for  computing  the  scale 
distance  parameter,  z. 


Case 

a 

b 

c 

1 

1 

4.3580 

1  7600 

0  0670 

2 

5.57"9 

1  7600 

0  0335 

3 

5.57~9 

1.7600 

0.0335 

A.3.2  Synoptic/METAR  Code.  Synoptic  total 
cloud-cover  reports  are  in  eighths.  METAR  does 
not  include  a  total  cloud-cover  report,  but  only 
inoividual  layer  amounts.  The  following  algorithm 
developed  at  USAFETAC  is  used  to  compute  the 
total  fractional  cloud-cover,  (T),  for  N  cloud- 
layers; 

T  =  Xy,  for  A/  =  1 

r=  0.85f;X^for  A/a2  08) 

M 


where  the  particular  value  of  Ay  is  determined 
from  Equation  (14)  (For  case  1,  Ay;,  etc.).  The 
coefficients  a,  b,  and  c  for  each  case  are  listed  in 
Table  A-1.  The  END  of  mean  cloud-cover  is 
computed  using; 


1.9962  +  0.(X)22e'^"' 

where  y  equals  O.Sfy^*,,  +  y,„)  for  Case  1  and 
simply  yi^  for  Cases  2  and  3.  The  sky-dome  scale 
distance,  (r),  is  then: 

r  =  ^  (17) 

2^ 

where  A  is  the  area  of  the  sky  dome,  usually 
assumed  to  be  2424  km^  (Burger,  1985). 


IF  T  -  1  THEN  T  =  1 

where  A',  is  the  fractional  cloud-cover  for  the  ith 
layer. 

The  relative  frequencies  of  each  of  the  categories 
(0/8,  1/8,  2/8,  ...,  8/8)  are  indicated  by  f,  where  i 
ranges  from  0  to  8  and  indicates  the  eighths  of 
coverage  for  that  frequency  interval.  The  ENDs  of 
/o  /<,  and  /,  coverage  are  computed  (y,^  and 
y;),  then  used  to  compute  the  following  set  of 
difference  values: 

(^y)i  =  /7  -  /o 

(^k)2  =  /.  -  /o  09) 

(A»3  =  /:  -  /4 

Two  additional  values  (p^  and  pjJ  are  computed: 
The  value  of  the  scale  distance  parameter  z  is 
computed  using  Equation  11  Equation  12  is  used 
to  compute  the  END  of  mean  cloud-cover,  {y^). 
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Pah  U  ~  U 
Pta-  fj  -  U 


(20) 


The  coefficients  used  in  this  equation  depend  on 
one  of  the  classes  (4,  5,  or  6).  In  most  instances. 
Case  4  applies.  Case  5  occurs  when  /,  <  and  /, 
<  and  Case  6  occurs  when/j  i  /,  and/j  < 

The  value  of  >  in  eq.  (4)  is  0.5(>7  +  for  Case  4 
and  for  Cases  5  and  6. 


TABLE  A*2.  Cloud-cover  distribution 
used  in  the  example. 


Rel 

Cum 

Cateaorv 

Freq 

Fteq  END 

CLR 

0.20 

0.20  -0.8418 

SCT 

0.25 

0.45  -0.1256 

BKN 

0.25 

0.70  0.5244 

OVC 

0.30 

1.00 

X.i.i  Example.  To  illustrate  how  these  variables 
are  computed,  consider  the  following  example. 
We  will  assume  a  location  that  uses  Airways 
reporting  has  the  following  distribution  of  cloud- 
cover  for  the  month  and  hour  of  interest:  clear 
20%,  scattered  25%,  broken  25%,  and  overcast 
30%.  The  cumulative  frequencies  and  their 
corresponding  ENDs  are  shown  in  Table  A-2. 


Xo  =  (0.5(0.5244  -  0.8414)1  ♦ 
(1.9962  ♦  (0.0022  x  exp(-0.89  x  4.9914)]} 
=  -0.0795 


(23) 


This  value  is  the  END  of  the  mean  sky  cover. 
This  corresponds  to  a  mean  sky-cover  of  0.53.  To 
compute  the  scale  distance: 


/-= 

2aMi4 


1.55 


(24) 


Thus,  the  Burger  Distribution  variables  for  the 
example  in  Table  10  are:  mean  sky  cover  0.53, 
and  scale  distance  1.S5  km. 

A.4  Calculation  of  the  Cumulative  Frequency 
Distribution 


A.4.1  Technique.  Next,  we  present  the  technique 
for  obtaining  the  continuous  cumulative 
distribution  of  sky-cover,  given  the  Burger 
distribution  variables.  To  determine  the 
probability,  (Pj,  of  obtaining  a  fractional  sky- 
cover  of  X  or  less,  we  compute  F  and  t: 


F  =  lOx.  X  s  0.5 
=  10(1  -  X),  X  >  0.5 


(25) 


/  =  1,  X  s  0.5 
=  -1,  X  >  0.5 


(26) 


This  example  follows  Case  1.  Computing  the 
initial  values; 

(Aji^  =  «« 

(0.5244  -  (-0.8418)1  =  1.3662 

Z  =  4.3580  ♦  (1.7600  x  In  (1.3362  ♦  0.0670)] 

^  4.8914 

(22) 


Recall  that  is  the  END  of  the  mean  sky-cover 
and  z  is  the  scale  distance  size.  We  then  carry  out 
the  following  calculations: 

y'o-ty,  (27) 

/V  =  5  -  INT(F )  (28) 

where  INT(F)  represents  the  integer  portion  of  F 
[e.g.,  INT(5.9)  “  5],  The  following  additional 
calculations  are  made: 
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=  (0.9981  ♦  0.0011  x  exp(0.89  x  r))  x  yj 

(28) 


The  END  of  the  desired  probability  is  -0.3104, 
which  corresponds  to  a  probability  of  38%.  Thus 
the  probability  of  obtaining  a  total  sky-cover  of 
30%  or  less  is  38%,  P{x  <  0.30)  =  0.38. 


IF  N  #  0; 

y,  =  >'c  *  [As  *  X  (C«  2)1 

X  (f  -  -  6) 

IF  N  >  1; 

w  (31) 

[At  -  B,  X  exp(C,  2)] 

/•1 

y,  =  fy.  (32) 

This  last  value,  (yj,  is  the  END  of  the  desired 
probability  P^. 


A.5  Calculations  with  the  Normal  Distribution. 
Since  the  integration  of  the  normal  probability 
density  function  has  no  known  analytical  solution, 
empirical  approximations  are  used.  These 
approximations  are  required  in  both  transforming 
into  and  out  of  the  Gaussian  domain.  The  latter 
calculation  is  straightforward. 

A.S.1  Converting  an  END  to  Its  Probability.  A 
fourth-order  polynomial  regression  is  used  to 
provide  an  approximate  solution  to  the  integration 
of  the  normal  probability  density  fuiiction.  To 
determine  the  probability  of  obtaining  an  END  of 
a  value  less  than  y,  ^(y),  the  following  algorithm 
is  used  (Abramowitz  and  Stegun,  1964): 


A.4.2  Example.  Continuing  with  the  example 
from  the  previous  section  (mean  of  0.53  and  scale 
distance  of  1.55  km),  what  is  the  probability  of  a 
sky-cover  of  30%  or  less  (x  “  0,30)? 

^  =  3  yo  =  -0.0785 
/  =  1  2  =  4.9914 

(33) 

yo  =  -0.0795 


IF  y  a  0  THEN  T  =  1 

(36) 

IF  y  <  0  THEN  T  =  -1 

fCy)  »  EH;{7y)/  (37) 

/.o 

The  coefficients  are  listed  in  Table  11. 


N  =  2 


y^  =  (1 .091 6)( -0.0795) 

=  -0.0868 

(34) 

Ya  =  -0.0868  ♦  (-1)  x 

(-0.0205  *  (0.0097  •  (0.52  •  4.9914)]) 

=  -0.1963 


y,  =  -0.1963  *0.0071  - 
(0.0086  X  exp(0.53  x  4.9914)]  (35) 

=  -0.3104 


TABLE  A-3.  Coefficients  used  to  convert 
from  a  probability  to  an  END. 


n 

-tL 

0 

1.189 

1 

0.2341 

2 

0.1370 

3 

0.0004091 

4 

0.023221 
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A.S.2  Converting  a  Probability  to  Its  END. 
Computing  the  END  (y),  from  a  probability  {P)  is 
more  involved.  The  technique  described  here  is 
from  Beasley  and  Springer  (1977).  Simpler,  but 
less  accurate  formulas  are  found  in  Boehm  (1976). 
First  we  compute  the  parameter  Q'. 

O  =  P  -  0.5  (38) 

If  IQ  I  s  0.42  use  Equation  39: 

y  =  -  (39) 

/•I 


Next  compute  the  parameter  S: 


S  =  [-ln(/7  )]’/* 


The  END  y  is  then: 


9-  T 


/•o 


1 


(41) 


(42) 


The  coefficients  y,  and  bj  are  listed  in  Table  A-4. 


The  coefficients  a,  and  Py  are  listed  in  Table  A-4. 
For  the  case  when  |Q|  >  0.42,  first  compute  the 
parameter  R\ 

IF  O  i  0  THEN: 

R  =  P 

(40) 

IF  O  <  0  THEN: 

/?  =1  -  P 


TABLE  A-4.  Coefficients  used  to  convert  from  a  probability  to  its  END. 


n 

Pn 

Yn 

0 

2.5066 

— 

-2.7872 

— 

1 

-18.6150 

-  8.4735 

-2.2980 

3.5430 

2 

41,3912 

23.0834 

4.8501 

1.6371 

3 

-25.4411 

-21.0622 

2.3212 

— 

4 

3.1308 
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ACRINABS 


BAA 

Burger  aerial  algorithm 

BKN 

broken  (cloud-cover  category) 

CFLOS 

cloud-free  line-of-sight 

CLR 

clear  (cloud-cover  category) 

CPCFLOS 

climatological  probability  of  cloud-free  line-of-sight 

OVC 

overcast  (cloud-cover  category) 

END 

equivalent  normal  deviate 

INT 

integer  function 

0-U 

Omstein-Uhlenbeck 

RTNEPH 

real  time  nephanalysis 

SCT 

scattered  (cloud-cover  category) 

SRI 

Stanford  Research  Institute 

SYT 

Simulation  and  Techniques  Branch 

USAFETAC 

United  States  Air  Force  Environmental  Technical  Applications  Center 

3DNEPH 

B-dime.tsional  nephanalysis 
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OD  4/OX.  Onixuka  AFB.  CA  94088-3430  .  1 

SSD  OO  4,  Onixuka  AFB.  CA  94088-3430  .  1 

Det  3.  Space  Systems.  Bldg  430.  Stop  77,  Buckley  ANQB.  CO  8001 1  -9599  .  1 

AFMC(I)/OOW.  Bldg  266.  Post  108P  Chidlaw  Rd..Wright-Patterson  AFB.  OH  45433-5000  .  1 

FASTC/TAW,  41 1 5  HebWe  Creek  Rd.,  Ste  33.  Wright-Patterson  AFB.  OH  45433-5637  .  1 

ASD/we,  Bldg  91 . 3fd  St  Wright-Patterson  AFB.  OH  45433-6503  .  1 

AFIT/CIR.  Wright-Patterson  AFB,  OH  45433-6583  .  1 

WL/DOA,  Wright-Patterson  AFB,  OH  45433-6543  .  1 

WiyOOW,  Wright  Patterson  AFB.  OH  45433-6543  .  1 

645  WS/CC.  Wright-pattersoo  AFB,  OH  45433-5000  .  1 

46  TQ/WE.  Holloman  AFB,  NM  88330-5000  .  1 

PUWE,  Kirtland  AFB,  NM  871 1 7-5000  .  1 

HQ  AFOTEC/WE,  Kirtlond  AFB,  NM  871 1 7-7001  .  1 

RL/WE.Qrlftiss  AFB,  NY  13441-5700  .  1 

RQMELAB/SUUCorridorW.  Ste  262.  26  Electronic  Pkwy.Qrlttiis  AFB,  NY  13441-4514  .  1 

AFCESA/WE,  Tyndall  AFB,  FL  32403-5000  .  1 

AFESC/RDXT.  Bldg  1 1 20.  Stop  21 .  Tyndall  AFB,  FL  32403-5000  .  1 

ESD/WE.  Vandenberg  Dr..  Bldg  1624.  Hanscom  AFB.  M A  01731-5000  .  1 

PL/TSML  Research  Library,  Hansccm  AFB,  MA  01731-5000  .  1 

PUGP.  Hansom  AFB,  MA  01 731  -5000  .  1 

3246TW/OOW,  Bldg  60,  Rm  60.  Eglin  AFB,  FL  32542-5000  .  1 

46  WS/CC,  EglinAFB.FL  31542-5000  .  1 

AFFTC/WE,  Edwards  AF3,  CA  93523-5000  .  1 

SMC/SDW.  PO  Box  92960.  Bldg  117,  El  Segundo,  Los  Angeles  AFB.  CA  90009-2960  .  1 

UTTR/WE.  Hill  AFB,  UT  64056-5000  .  1 

USCENTCOM/CCJ3-W,  Bldg  540,  MacDIII  Blvd,  MacDill  AFB.  FL  33608-7001  .  1 

HQ  AFTAC/TNK,  1030  Highway  A1  A,  Patrick  AFB.  FL  32925-5000  .  1 

ESMC/WE,  Patrick  AFB,  FL  32925-5000  .  1 

OL-A.  AFCOS,  Site  R,  Fort  Ritchie.  MO  21719-5010 .  1 

USAFALCENT  RA.  Pope  AFB,  NC  28308-5000  .  1 

CCSO/FL,  Tinker  AFB,  OK  73145-6340  .  1 

AFOSR/NU  Bolling  AFB.  DC  20332-5000  .  1 

TFWC/WE,  Nellis  AFB,  NV  89191-5000  .  1 

SMC,  Det  2/TDO,  Onlluke  AFB  CA  94088-3430  .  1 

AL/OEBE.  2402  East  Drive,  Brooks  AFB,  TX  78235-51 14  .  1 

AMC/XOWR,  Bldg  P40  N,  Martin  Ave.  Soott  AFB.  IL  62225-5000  .  1 


23 


HQ  AFSOC/DOW,  Bldg  1,  Hurtburt  FLO,  FU  32544-5000 .  1 

ATC/DOTY/,  244  F  Str**l  East,  Suita  #3,  Randolph  AFB,  TX  78150-4325  .  1 

Oat  5,  HQ  AWS.  Kaaslar  AFB.  MS  39534-5000  .  1 

CFA,  C-2/SWO.  APO  AP  96258-0210  .  1 

PACAF/DOW,  Bldg  1102,  Hickain  AFB.  HI  96853-5000 .  1 

Oat  1 .  HQ  PACAF.  COMNAVMAR,  PSC  489.  Box  20.  FPO  AP  96540-0051  .  1 

11WS,  1215  FlighBina  Ave.  Ste  2,  Eielson  AFB.  AK  99702-1520  .  1 

USSTRATCOM/J3615,  Rm  LI  27,  Bldg  522,  901  SAC  Blvd.  Offutt  AFB.  .ME  68113-5000  .  1 

2  WS/CC.  Bldg  5546.  245  Davis  Ava.  Barksdala  AFB,  LA  71 1 1 0-5002  .  1 

ACC/DOW.  Bldg  21 .  30  Elm  St  Sta  21 5.  Langlay  AFB.  VA  23655-2093  .  1 

24WS,  Unit  0640,  APO  AA  34001-5000  .  1 

9COS/AOSW.  Bldg  1130,  Shaw  Df.,  Shaw  AFB.  SC  29152-5410  .  1 

ATSI/CDW,  US  Army  Intal.  Ft  Huachuca  Al.  AZ  85613-5000  .  1 

USCENTCOM/CCJ3-W,  MacOill  AFB.  FL  33608-7001  .  1 

USSOCCENT/SCJ2-SWO.  MacOlll  AFB.  FL  33608-7001 .  1 

USSCX:OM/SOJ3-W,  MacOill  AFB.  FL  33608-6001  .  1 

1WG.  Bldg  168,  Hardae  St.  Ft  McPharson.  GA  3030C'-5000  .  1 

USAFE/DOW,  Unit  3050,  Box  15.  APO  AE  09094-5000  .  1 

17AF/DOW.  Unit  4065,  APO  AE  09136-5000  .  1 

HQ  USEUCOM  ECJ3.  Unit  30400.  Box  1000,  APO  AE  09128-4209  .  l 

7WS,  CINCUSAREUR/AREAWX.  APO  AE  09403-5000  .  1 

16AF/SWO.  Unit  6365.  APO  AE  09641-5000  .  1 

COMNAVOCEANCOM.  Coda  N312.  Stannis  Spaca  Ctr.  MS  39529-5000  .  2 

COMNAVOCEANCOM  (Capt  Brown,  Coda  N332),  Stannis  Space  Ctr.  MS  39529-5001  .  1 

NAVOCEANO  (Rusty  Russum),  Bldg  8100,  Rm  2030.  Stannis  Spaca  Ctr.  MS  39522-5001  .  2 

NAVOCEANO,  Coda  9220  (Tony  Oftolano).  Stannis  Space  Ctr.  MS  39529-5001  .  1 

Maury  Oceanographic  Library.  Naval  Ooaarwgraphy  Offioa,  Stannis  Spaca  Ctr,  MS  39522-5001  .  1 

Naval  Research  Laboratory,  Monterey,  CA  93943-5(X)6  . .  1 

Naval  Research  Laboratory,  Coda  4323,  Washington,  DC  20375  .  1 

Naval  Research  Laboratory  (Dr  Riley),  Coda  4180,  Washington,  DC  20375  (SESS  reporte) .  1 

Naval  Postgraduate  School,  Chmn,  Dept  of  Meteorology,  Coda  63.  Monterey,  CA  93943-5000  .  1 

Commanding  Olfioar,  Naval  Polar  Oceanography  Canter.  4301  Soitland  Road,  FOB  #4,  Washington,  DC  20395-5108  .  1 

Naval  Oceanography  Command  Ctr.  COMNAVMAR  Box  12,  FPO  San  Francisco,  CA  96630-5000  .  1 

Commanding  Olfioar.  Naval  Oceanography  Command  Ctr.  PSC  819,  Box  13.  FPO  AE,  09645-3200  .  1 

Chief,  APG  Met  Team,  Bldg  1134.  Attn;  AMSTE-TC-AM  CAB,  Aberdeen  Proving  Ground.  MD  21 005-5001  .  1 

Atmospheric  Sdencas  Laboratory  (SLCAS-AS-I  3  10-2c),  White  Sands  Missile  Range,  NM  88002-5501  .  1 

TECOM  Atmos  Sd  Div,  AMSTE-TC-AA  (MacBlain),  White  Sands  Missile  Range,  NM  88002-5504  .  1 

White  Sands  Met  Team,  AMSTE-TC-AM  (WS),  White  Sands  Missile  Range,  NM  86002-5501  .  1 

USA  TECOM,  ATTN:  AMSTE-TC-AM  (RE)  TECOM  Met  Team,  Redstone  Arsenal,  AL  35898-8052  .  1 

Director,  U.S.A.-CETEC.  Attn:  GL-AE  (Whitmarsh),  Fort  Belvoir.  VA  22060-5546  .  1 

Technical  Ubrary,  Dugwsy  Proving  Ground,  Dugway,  UT  84022-5000  .  1 

HQ  NATO  Staff  Metaorol^ical  Officer  IMS/OPS  APO  AE  09724  .  1 

NOAA/MASC  Ubrary  MC5,  325  Broadway.  Boulder.  CO  80303-3328  .  2 

NOAA  Ubrary.EOC4W5C4.  Attn;  ACQ,  6009  Executive  Blvd.  Rockville.  MO  20852  .  1 

NOAA/NESDIS  (Attn;  Nancy  Everson,  E;RA22),  World  Weather  Bldg.  Rm  703,  Washington,  DC  20233  .  I 

NGDC,  NOAA,  Mail  Code  E/GC4,  325  Broadway,  Boulder.  CO  80333-3328  .  1 

NWS  W/OSD,  Bldg  SSM  0-2  East-West  Hwy,  Silver  Spring,  MO  20910  .  1 

NIST  Puba  Production,  Rm  A835,  Admin  Bldg,  Gaithersburg,  MD  20899  .  1 

NASA-MSFC-ES44,  Attn;  Dale  Johnson.  Huntsville,  AL  35812-5000  .  1 

NASA-MSFC-ES44,  Attn;  Qwenevere  Jasper,  Huntsville,  AL  35812-50(X> .  1 

OTIC-FDAC.  Canneron  Stobon,  Alexandria,  VA  22304-6145  .  2 

AUL/LSE,  Maxwell  AFB,  AL  361 12-5564  .  . .  1 

AWSTL  Scotl  AFB,  IL  62225-5438  .  35 
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